NMR relaxation time around a vortex in stripe superconductors 
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Site-dependent NMR relaxation time T\ (r) is calculated in the vortex state using the Bogoliubov- 
de Gennes theory, taking account of possible "field-induced stripe" states in which the magnetism 
arises locally around a vortex core in d-wave superconductivity. The recently observed huge en- 
hancement T 1 ~ 1 (r) below T c at a core site in Tl2Ba2CuOe is explained. The field-induced stripe 
picture explains consistently other relevant STM and neutron experiments. 

PACS numbers: 76.60.Pc, 74.60.Ec, 74.25.Jb 



Much attention has been focused on microscopic elec- 
tronic structure, or low-lying excitation spectrum around 
a vortex core in a superconductor, particularly in the 
high T c cuprate superconductor, whose pairing mecha- 
nism is still not known. It is expected that the vortex 
core might be a useful window through which we can 
probe the ground state in the " normal phase" by remov- 
ing superconductivity. 

Recently, keen interest is stimulated by a series of re- 
markable experiments: Neutron scattering experiments 
observe an enhancement of the incommensurate mag- 
netic satellite peaks by applying magnetic fields [jl], ||, 
|j. The checkerboard pattern of the local density of 
states (LDOS) in E^S^CaC^Og+a, or the so-called 
"halo" is found around a vortex site by scanning tun- 
neling microscopy (STM) experiments by Hoffman et 
aZ.Q]. Vortex-imaging NMR experiments detect anoma- 
lous relaxation time T\ around a vortex core site in 
YBa 2 Cu 3 07@ and YBa 2 Cu 4 8 §]. The set of these ex- 
periments collectively point to a notion that incipient 
spin density wave and/or charge density wave, namely, 
stripe states in a broad sense are induced around a vortex 
core by applying a moderate magnetic field. 

In d-wave pairing, which is now fairly well established 
for high T c cuprate superconductors, the low-lying vortex 
core excitations are expected to have a finite density of 
states (DOS) exactly at the Fermi level. Thus the core 
should be filled up by these quasi-particles and tunnel- 
ing conductance as a function of the bias voltage should 
exhibit a zero-energy peak (ZEP) structure at a core. 
However, the existing STM results Q] so far fail to detect 
this ZEP, instead they see the empty core, namely there 
is not any clear structure in the tunneling conductance 
at a core site. This disparity between the theoretical pre- 
diction and experimental facts is rather serious, because 
in a conventional superconductor such as 2H-NbSe2 the 
observed quasi-particle spectra by Hess et aZ.f| are al- 
most perfectly understood theoretically in terms of the 
bound state trajectories ||. The situation becomes more 
serious since the vortex imaging NMR experiments done 
by Mitrovic, et a/.|| and Kakuyanagi, et aZ.^J found that 
the site-dependent T\ becomes longer at the vortex core, 
after being shorter with approaching towards the vor- 
tex. It is contradicted with existence of the ZEP of the 
LDOS, which should lead to the shortest T\ at the core. 



This suggests absence of the ZEP associated with d-wave 
vortex, being rather consistent with the STM0. 

The emerging notion of the field induced stripe 
matches with the above neutron experiments because the 
four site charge modulation corresponding to the STM 
observation^) is accompanied by the eight site spin den- 
sity wave, or stripe, seen as the incommensurate satel- 
lites. The applied field with a mere few Tesla enhances 
those satellite intensities and the estimated induced mo- 
ment is roughly proportional to -ffjjj. 

Prior to the present work, there have been several the- 
oretical attempts to explain some of these facts, notably 
the STM experiments under |l| |l| g |l§, 

PI , |l7[] and without i?Jl|, These theories are based 
on the field induced incipient spin and/or charge density 
orders, giving an explanation of the STM observation of 
four-site periodic LDOS modulation. 

Quite recently, two new NMR imaging experiments |2C|, 
pl| on the T-dependence of T\ are reported: Kakuyanagi, 
et oZ.@ show in Tl 2 Ba 2 Cu0 6 that as lowering T, l/T x T 
at a core site increases divergingly towards a tempera- 
ture Tm below T c and then decreases at lower T . This 
is contrasted with \jT{T at other sites which exhibit a 
monotonic decrease below T c . This result suggests that 
below Tm the magnetism appears locally exclusively at a 
vortex core site and other sites stay in the normally ex- 
pected d-wave state. Mitrovic, et aZ.[^l| observe a similar- 
divergent behavior in \jT\T of YBa2Cu307 for the core 
site where the crossover temperature Tm is quite low and 
\jT{T for the sites outside the core shows a constant at 
low T . The implication of these experiments is two-fold: 
The missing ZEP in d-wave vortex core must be present 
in the limited T and H region. Below Tm the locally field 
induced stripe must exist, which is to remove the ZEP, 
giving rise to a suppression of T{~ at lower T. 

The purpose of this study is to demonstrate that 
the field induced stripe picture explains these experi- 
mental facts in a coherent way. This is done by ex- 
plicitly calculating the site-dependent T±, based on the 
Bogoliubov-de Gennes (BdG) equation for d-wave super- 
conductors. Our calculation shows that the vortex imag- 
ing NMR(2^, Q yields a wealth of information on 
d-wave superconductors under H. 

We begin with the standard Hubbard model on a two- 
dimensional square lattice, and introduce the mean field 
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ni,o = ( a la a i,<?) at tne *-site, where a is a spin index and 
i = (i x ,i y )- We assume a pairing interaction V between 
nearest-neighbor (NN) sites. This type of pairing inter- 
action gives d-wave superconductivity [^2| |23| ? ]. Thus, 
the mean-field Hamiltonian under H is given by 

ft = - hj a i* a j,v + U Y1 ni- a a\ >cr a ita 

i,j.(T i.<r 
+ V J2 ( A L, CT «*,- a *+e- + ^,i,A,A+e,-.) (!) 

where a\ (ai j(T ) is a creation (annihilation) operator, and 
i + e represents the NN site (e 
integral is expressed as 



±i, ±y). The transfer 
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FIG. 1: Spatial structures of the d-wave order parameter 
A,j(r) and the induced staggered moment M(r). At the vor- 
tex core site V, |A d (r)j = 0. The unit cell size is 24 x 24. 
U/t = 3.0. T/t = 0.01. 



ti,j = k,j ex P[i^- / A ( r ) • dr ]j 



(2) 



with the vector potential A(r) 
gauge, and the flux quantum i 



x r in the symmetric 
For the NN pairs 
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)~ AijSj.i+e and a is an index of the eigenstate|22[ 
The self-consistent condition for the pair po- 
tential and the number density is given by A,j = 

( a jd a iA) = E a «a( r «K( r j)/( £ a). n i,1 = ( a !,T ai <T) = 

E a \u a (ri)\ 2 f(E a ), n u = (at ; a u ) = J2 a K(n)\ 2 {l - 
f(E a )). The charge density n, = nj^+n^j., the spin den- 
sity S z ,i = —1H,\) and the staggered magnetization 
Mi = S z i(— \) l x +l y . The d-wave order parameter at site 
% is A d>i = (A iti + A_z }i - Ay^ - A^y :i )/4 with 
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A; 



i+ i exp[i- 



(r;+r 



i)/2 



A(r) ■ dr], 



where A i<i+S = (aj+gja^f) - (a l+i ^a hl ). 

We typically consider the case of a unit cell with 24 x 24 
sites, where two vortices are accommodated. The spa- 
tially averaged hole density is set to = 1 — nl ~ | 
by tuning the chemical potential /i. By introducing the 
quasimomentum of the magnetic Bloch state, we obtain 
the wave function under the periodic boundary condition 
whose region covers many unit cells. 



We construct Green's functions from E a , u a (r), 
v a (r), and calculate the spin-spin correlation function 
X+,- ( r ; r 'i i&n) p3|] ■ We obtain the nuclear spin relax- 
ation rate. 



Uj — t. For the next-NN pairs situated on a diagonal R(r,r') = 
position on the square lattice, t<j = t'. For the third-NN 
pairs, which are situated along the NN bond direction, ~~ 
tij = t". To reproduce the Fermi surface topology of 
cuprates, we set t' = -0.12i and t" = 0.08*||f|. We 
consider the pairing interaction V = 2.0t. The essential 
results of this paper do not significantly depend on the 
choice of these parameter values. 

In terms of the eigen-energy E a and the wave functions 
u a (ri), v a {vi) at the i-site, the BdG equation is given by 



Imx+,_(r, r', tfl n - SI + i V )/(n/T)\ n ^ 
-5> a (rK(rOMrK,(r') 



-u a (r)<(r> a( (r)«*,(r')] 
xirTf{E a )5(E a -E a ,) 



(4) 



with the Fermi distribution function f(E). We consider 
the case r = r' by assuming that the nuclear relaxation 
occurs at a local site. Then, r-dependent relaxation time 
is given by Ti(r) = l/i?(r,r). We use S(x) = 7r _1 Im(a: — 
irj)^ 1 to handle the discrete energy level of the finite size 
calculation. We typically use rj = O.OTi. In Eq. (0), 
the first term is proportional to the N^(r, E) x Ni(r, E) 
when r = r'. To understand the behavior of 2~i(r), we 
also consider the LDOS given by N(r,E) = N^(r,E) + 
Ni(r, E) = £JK(r)| 2 <5(£ Q - E) + \v a (r)\ 2 S(E Q + E)]. 

By solving Eq. (||) self-consistently, we obtain the order 
parameter A^i and the sublattice magnetization Mi as 
shown in Fig.l. It is seen that the moment is induced 
exclusively around the vortex core where the order pa- 
rameter vanishes. The checkerboard modulation of M 
with eight site period along x and y directions is su- 
perimposed. This is barely seen in Fig. 1(b) where the 
"floor" in M is modulated with eight-site period. In 
the Fourier transformation of S z ^ the eight-site period 
spin structure gives rise to peaks at the corresponding 
wave number Q in addition to the peak of the vortex 
lattice period. As U increases the checkerboard modula- 
tion with two ordering vectors Q = (1^,0) and (0, §7r) 
changes into a one-dimensional stripe form character- 
ized by Q = (§7r,0) or (0, §tt) (see Fig.2 in Ref.[@ for 
stronger stripes). The induced moment at the core in- 
creases also with U. This spin modulation accompanies 
the charge modulation with 2Q, namely, four site peri- 
odicity. The core accommodates excess particles, which 
is opposite to the stripe-free case where it does excess 
holes [^6|. These results basically coincide with the previ- 
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FIG. 2: Local density of states iV(r, E) at the core site for 
U/t = 0, 2.6, 3.0. The zero-energy state at E = for U/t = 
is gradually suppressed with U. Inset shows the spatial profile 
of the zero-energy LDOS N(r, E = 0) for U/t = 3.0. T/t = 
0.01. 

ous studies JlJ, [01 [l4|, [l^, [ll| and also capture the essen- 
tial features observed by STM experiment Q . The vortex 
core changes to a circular shape by induced moment as 
seen from Fig.l (a). When U = its shape was square 
like. H 

In Fig. 2 we plot the LDOS as a function of the en- 
ergy E/t at the core site. For the small U/t = 2.6 case 
where the appreciable moment is not induced, the ZEP 
appears, centered at E/t — 0. As U/t increases, the ZEP 
splits and is gradually suppressed. A gap-like feature be- 
comes evident for U/t — 3.0 in Fig. 2. It is noted that the 
removed DOS at E — piles up inside the bulk gap at 
\E/t\ ~ 1. This DOS feature resembles the bound state 
at E =8meV observed by STM on YBa 2 Cu 3 7 §. The 
inset shows the spatial distribution of the zero-energy 
state: When the induced moment is absent (U = 0), 
the zero-energy state accumulates at the core site. This 
is altered by introducing the induced moment, namely, 
iV(r, E = 0) decreases towards a core after taking a max- 
imum. This is contrasted with the peak structure at a 
core in the stripe-free case (see Fig. 6 in Ref . [p3| ) . We also 
note that the LDOS at E — shown in inset contains the 
Fourier component with 4-site checkerboard modulation. 

We depict the spatial distribution for T 1 _1 (r) at low 
T(= 0.01*) for U/t = 3.0 in upper panel of Fig.3. It is 
seen that corresponding to the depressed LDOS at the 
core shown in Fig. 2 due to the induced moment, T 1 ~ 1 (r) 
diminishes around the vortex core after it increases with 
approaching the vortex from far sites. In lower panel 
of Fig.3, T 1 _1 (r) is depicted along the selected paths in 
the vortex lattice shown in upper panel. It is clearly seen 
that (1) T 1 _1 (r) at the core site V takes a local minimum. 
(2) One of the minima of T 1 _1 (r) occurs at C-point. (3) 
The relative T 1 _1 (r) values at V and C points depend 
on the magnitude of the induced moment. When \M\ is 
large enough, T-j -1 ^) becomes the absolute minimum at 




FIG.3: Spatial profile of Tj~ (r) (upper panel). T 1 _1 (r) varia- 
tion along the path (lower panel): V^C^S-^V where empty 
(filled) squares correspond to U/t — 3.0 (U/t — 0). In the 
stripe- free case (U/t = 0) Tf^r) at V is extremely large. T? 
is the normal state value at T c . 

V. (4) The maximum T 1 ~ 1 (r) occurs near the core site 
whose distance from the core depends on \M\. In larger 
\M\ case this distance becomes large because of the larger 
LDOS suppression due to the induced moment. 

Both site-selective experiments of 17 NMR on 
YBa 2 Cu 3 7 § and YBa 2 Cu 4 8 § for the Cu0 2 planar 
site show that T 1 _1 (r) at the core is smaller than the 
neighboring sites for a certain H and T region. The 
Tj~ (r) maximum occurs in the vicinity of the core, i.e. 
T 1 _1 (r) decrease with increasing the internal field around 
the vortex core region. These facts are reproduced by our 
calculation, evidencing the usefulness of the site-selective 
NMR method in general. When we see the Tf^r vari- 
ation along C-^S^V in the experimental data]|[ [ij. 
T± (r) has a minimum near the saddle point S. This be- 
havior is also reproduced in our calculation of the stripe 
state here, while in the previous calculation (U = 0), 
T 1 _1 (r) monotonically increases along C— >S-^V as dis- 
played in lower panel in Fig.3. 

We show the T-dependence of l/Ti(r)T for three sites 
in Fig. 4. It is seen that l/7i(r)T at the vortex site V 
displays a large enhancement around Tm below which 
M begins to appear as shown in the inset. Thus for 
Tm <T <T c the zero-energy state increases 1/T\T and 
for T < Tm the induced moment opens a gap at E = 
as shown in Fig. 2, reducing 1/T{T towards low T. The 
\/T{T behavior in the nearest neighbor site to the core is 
moderately enhanced at Tm while that in S-site remains 
same as in the stripe-free case. These correspond to the 
fact that the induced moment is limited only near the 
core. We should remark here that the enhancement in 
\/T{T below T c delicately depends on applied field. As 
H increases, it is suppressed rather easily (see Fig. 12 in 
Ref.@). 

The diverging behavior in l/Ti(r)T at the core site 
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FIG. 4: T-dependence of l/Ti(r)T for three sites. (1) 
The vortex core (V). (2) The nearest site to the core. (3) 
The saddle point site (S). Inset shows the corresponding In- 
dependence of Ad and M at their maximum values. U/t — 
3.0. 

is strikingly similar to that observed in 205 Tl NMR on 
Tl 2 Ba 2 CuO6^0|. Our calculation is also consistent with 
1/Xi (r)T behavior at the S point: Their data show strong 
decrease of l/Ti(r)T at S below T c and l/TiT=constant 
at lower T, indicative of the zero-energy state. 

As for Mitrovic's data[fl| , they observe that (A) 1 /T{T 
inside the core is enhanced as lowering T while (B) 1/T\T 
outside the core stays at a constant which increases as 
H increases. These behaviors are understandable if we 
assume that Tm is low in YBa2Cu307, namely their ex- 
periments are done for T > Tm : We can explain the fact 
(A) because as lowering T the thermal smearing factor 
coming from f'(E) in eq.(^) becomes narrow enough to 
see the sharp ZEP, driving \jT{T to increase towards 
T — > 0. The fact (B) is consistent with the previous 
more detailed NMR at the saddle point by Zheng, et 



al.fZ^ on YBa 2 Cu 4 8 : They conclude that 1/TiT oc H 
at the saddle point at low T. As shown previously the 
spatial average DOS N(0) oc H 0A1 for d-wave casep8|, 
then 1/TiT oc iV(0) 2 - H S2 , explaining the approxi- 
mate TJ-linear behavior within the present experimental 
accuracy. 

Many ordinary T± measurements are done by using the 
resonance field at the maximum intensity which corre- 
sponds to the saddle point S of the field distribution in 
the vortex state. Above T c , 1/TiT shows a variety of 
different T-behaviors, depending on probed nuclei, such 
as O, Cu, Y, etc, or materials, or doping levels. These 
concern the so-called pseudo-gap phenomena. In con- 
trast, the low T behaviors, remarkably enough, are quite 
similar, namely, a T 3 -law below T c which is crossed over 
to the T-linear law at lowest T. The former T 3 -law in- 
dicates a line node for d-wave pair while the T-linear 
law is indicative of the zero-energy state outside of the 
vortex core, in spite of possible existence of the induced 
stripe. This is due to the fact that the induced moment 
is confined in the narrow core region and LDOS or T^ 1 is 
modified only there as seen from the lower panel of Fig. 3. 

In conclusion, we have calculated the site-dependent 
relaxation time T 1 _1 (r) for NMR, taking into account 
the possibility of the field- induced moment. The large 
Tj~ (r) enhancement observed in the vortex core site 
NMR in 205 T1 of Tl 2 Ba 2 Cu0 6 is explained by our calcu- 
lation. Other experiments for the vortex imaging^], ^, |2l[] 
are also analyzed in a coherent manner. An emerg- 
ing notion through these analyses, together with other 
STM and neutron scattering experiments, is the field- 
induced stripe picture in the vortex state of high T c 
superconductors. The strength of the induced stripe 
depends on material's parameters. The existing data 
show that Tl 2 Ba 2 CuOg is strong while YBa 2 Cu307 and 
YBa 2 Cu40s are weak. 

We thank K. Kumagai, Y. Matsuda and K. Ishida for 
useful discussions on their NMR. 
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